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ABSTRACT 


This  first  quarter  the  reflection  and  scattering  models  were 
further  developed,  and  a review  of  the  literature  preliminary 
to  the  development  of  a convergence  zone  propagation  model 
was  completed. 

In  Section  II,  the  expansion  of  the  heuristic  scattering  model 
to  include  composite  (both  large  and  small  scale  roughness) 
surface  statistics  is  given  which  should  improve  backscattering 
predictive  capabilities.  This  model  does  not  include  slope 
correction  terms  but  has  the  advantage  of  simplicity.  The 
more  rigorous  theoretical  model  which  includes  slope  correc- 
tions terms  is  also  expanded  to  allow  for  composite  surface 
roiaghness . 

A survey  of  the  literature  covering  convergence  zone  propagation 
is  summarized  in  Section  III.  An  approach  is  outlined  for  the 
development  of  a convergence  zone  model  using  a Green's 
function  expansion.  This  technique  results  in  an  ordered 
integral  expansion  where  each  integral  can  be  associated  with 
a particular  multipath. 
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I . INTRODUCTION 

Work  done  under  Contract  N00024-71-C-1266  for  the  period  of 
1 April  1972  to  30  June  1972  is  simmarized  in  this  report.  This  vork 
represents  a continuing  investigation  of  the  reflection  and  scattering 
of  sound  from  the  rough  ocean  boundaries,  and  the  interaction  of 
scattering  with  sound  propagation  in  the  ocean.  In  Section  II  the  two 
scattering  models  derived  at  ARL  are  extended  to  apply  to  composite 
surfaces.  This  work  will  provide  a basis  for  the  prediction  of  back- 
scattering  from  the  ocean  surface.  It  will  be  shown  that  both  large 
and  small  scale  surface  roughness  effects  are  important  in  the  back- 
scatter  case.  In  addition,  work  on  convergence  zone  propagation  was 
begun  this  quarter.  A brief  survey  of  the  literature  was  completed 
and  is  given  in  Section  III.  Theoretical  treatments  of  convergence 
zone  phenomena  are  discussed  for  both  ray  and  wave  theories . The 
limitations  of  both  classical  ray  and  wave  methods  are  summaxized. 

It  was  found  that  limitations  of  normal  mode  theory  can  be  overcome 
by  transforming  the  basic  contour  integral  into  an  infinite  sum  of 
multipath  integrals.  The  evaluation  of  these  ray  integrals  is  briefly 
discussed  with  an  outline  of  future  work  dealing  with  this  approach. 
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II.  BACKSCATTERING  FRCM  COMPOSITE  SURFACES 
A.  Introduction 

During  this  quarter  the  two  scattering  models  developed  at  ARL 
were  revised  and  extended  to  give  theoretical  predictions  of  the 
backs cattering  from  composite  rough  surfaces.  The  simultaneous 
development  of  both  scattering  models  was  carried  out  in  order  to 
better  understand  the  importance  of  phase,  beam,  and  slope  approxi- 
mations in  the  backscatter  case.  It  was  found  that  the  exact  forms 
of  the  phase  and  beam  approximations  are  relatively  unimportant  as  long 
as  the  coherence  of  the  scattered  field  near  normal  incidence  is 
properly  accounted  for.  However,  it  was  also  found  that,  for  the 
backscattering  case,  the  exact  slope  treatment  must  be  used,  and  that 
the  effects  of  shadowing  should  also  be  included. 

The  scattering  coefficients,  a,  used  in  this  section  are  related 
to  the  usual  definition  of  scattering  strength  through  S=^a/A,  where 
A is  the  active  scattering  area.  Here  it  is  assumed  that  the  scattering 
strength,  S,  is  referenced  to  1 yd.  In  the  following  section  the 
heuristic  (stationary  phase)  scattering  model  is  developed  from 
separate  calculations  of  the  coherent  and  incoherent  parts  of  the 
scattered  field.  The  modified  Fresnel  model  with  the  exact  slope 
treatment  is  given  in  Section  C and  is  compeired  with  the  heuristic 
model. 


B. 


Heuristic  Model 


The  heuristic  scattering  model  (given  in  detail  in  the  final 
report  under  Contract  N0002U-69-C-1275  aJid  in  "Forward  and  Specular 
Scattering  from  a Rough  Surface;  Theory  and  Experiment",  by  Boyd  and 
Deavenport^)  is  here  revised  and  extended  to  apply  to  composite  rough 
surfaces.  Particular  emphasis  is  placed  on  improving  the  backscatter- 
ing  predictions  without  degrading  the  forward  and  specular  scattering 
predictions.  This  subdivision  will  be  divided  into  three  parts 
following  the  chronological  development  of  the  model.  In  the  first 
part  a detailed  derivation  of  the  coherent  pressure  for  composite 
surfaces  will  be  given.  The  second  part  will  give  a brief  discussion 
of  the  composite  characteristic  functions  used  in  the  scattering 
model.  The  third  part  will  contain  the  development  of  the  incoherent 
scattering  coefficients  for  various  surface  statistics  and  some 
numerical  results  for  the  backscatter  case . 


1.  Coherent  Pressure 

The  scattered  pressure  is  given  by 

e^^  . e^dxdy  .( l) 

The  various  constants  and  ranges  appearing  in  Eq.  (l)  are  defined 
in  Ref.  1.  Assvme  that  the  surface  height  is  given  by 
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where  ^^(x,y)  and  ^^(x.y)  are  random  variables  representing  the  laxge 
and  small  scale  components  of  the  siarface  roughness,  respectively. 

Eq . ( 1 ) then  becomes 
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The  usual  development  of  the  coherent  pressure  can  now  be  made  using 
this  equation.  The  mean  pressue  is  given  by 
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, etc.,  and 


is  the  six-dimensional  probability  density  of  surface  heights  and 
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slopes.  Substituting  Eq.  (3)  into  Eq.  (4),  interchanging  the  order 
of  integrations,  and  using  the  definition  of  the  characteristic 
function  gives 
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If  the  slopes  heights  assumed 

independent , then 
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Further,  if  the  slopes  ^]^x’^2x’^ly’^2y  zero-mean  processes, 

then  the  expectation  value  becomes 
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Thus,  if  ^^(x,y)  and  ^2(x,y)  are  isotropic,  the  remaining  expectation 
value  can  be  removed  from  the  integral  and 
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where  the  term  in  braces  is  the  pressure  reflected  from  a plane 
surface.  The  coherent  intensity  is  then  given  by 


^s^COH 
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where  is  the  plane  surface  intensity  and  the  other  factor  is  the 
square  of  the  composite  characteristic  function  associated  with  the 
distribution  of  surface  heights.  The  coherent  scattering  coefficient 
is  obtained  from  this  by  dividing  out  the  intensity  scattered  from 
a plane  surface  in  the  specular  direction. 
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2.  Composite  Characteristic  Functions 


The  one-dimensional  composite  characteristic  function  used 
in  the  scattering  model  is  given  by 

/-ikr(L+Co)\ 


If  and  ^2  independent,  this  becomes 
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where,  for  a Gaussiein  process 
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and  for  an  exponential  process 
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Here  g=k  y h , where  h is  either  the  large  or  small  scale  rms  height. 


The  two-dimensional  composite  characteristic  functions  can 
be  treated  similarly  if  the  large  and  small  scale  processes  are 
independent : 

/ik7[(^^-H^2)  - /ik7(q-qx 


For  the  Gaussian  bivariate  distribution 
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and  for  the  exponential 
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5.  Incoherent  Scattering  Coefficient 


The  incoherent  scattering  coefficient  is  given  by 
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An  initial  attempt  at  adapting  this  expression  to  composite  surfaces 
consists  of  simply  replacing  the  characteristic  functions  with  their 
equivalent  composite  characteristic  functions,  i.e., 


(^8) 


It  should  be  noted  at  the  outset,  however,  that  this  approach  is  an 
approximation  in  that  no  modification  is  made  to  the  slope  terms 
which  appear  in  the  original  integral  for  the  total  pressure.  These 
slope  terms  have  already  been  evaluated  at  the  stationary  phase 
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points  and  removed  from  the  inte^/ral  as  Beckmc-.in's  F.  The  effect  of 
this  approximation  will  be  examined  in  the  following  section. 

Equation  (17)  will  here  be  developed  for  (a)  Gaussian  large 
and  small  scale  distribution  functions , (b)  exponential  large  and 
small  scale  distributions,  and  (c)  mixed  Gaussian  and  exponential 
distributions.  Forms  for  both  Gaussian  and  exponential  correlation 
functions  will  be  given  in  all  three  cases. 


(a)  For  Gaussian  large  and  small  scale  processes,  Eqs . (ll).  (l2), 
(l4),  (15),  and  (18)  are  used  so  that  Eq.  (17)  becomes 
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and  h^,  h^  are  the  large  and  small  scale  rms  heights.  To  evaluate 


Eq.  (19)  the  characteristic  function  must  be  expanded  in  terms  of  the 
Jarpe  and  small  scale  correlation  functions  and  Cg: 
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where  the  firot  term  is  required  to  be  unity: 
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where  the  J=0,  n=0  term  is  1.  Thus  the  integral  becomes 
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where  the  o=0>  '^=0  term  of  the  sum  is  omitted.  The  integral  to  be 


evaluated  i 


where  and  are  the  large  and  small  scale  correlation  functions, 
respectively.  For  Gaussian  correlation  functions 
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where  and  are  the  large  and  small  scale  correlation  lengths. 
Equation  (2J)  "becomes 
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Performing  this  integral  and  substituting  the  result  into  Eq.  (22) 
gives  (for  Gaussian  distributions  and  correlation) 
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For  exponential  correlation  functions 
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so  that  Eq.  (2j)  becomes 
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Thus  for  Gaussian  large  and  small  scale  iistribution  functions  and 
exponential  large  and  small  scale  correlation  functions,  the  inco- 
herent scattering  coefficient  is 
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(b)  For  exponential  large  and  small  scale  distribution  functions 
Eqs.  (ll),  (15),  (1^),  (16),  and  (18)  are  used  in  Eq.  (l?)  so  that 
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The  characteristic  f-unction  expansion  is  given  by 
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Again,  the  first  term  of  the  expansion  is  equal  to  the  one-dimensional 
characteristic  function  so  that  Eq.  (30)  becomes 
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where  the  n=0,  j=0  term  is  omitted. 


13 


M 


For  Gaussian  correlation  functions  Eq.  (25)  is  used  so  that 
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For  exponential  correlation  functions  Eq.  (28)  is  used  in  Eq.  (J2)  and 
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(c)  Aspuming  independent  large  and  small  scale  processes,  the 
integral  in  Eq,  (l?)  becomes 
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Assuming  that  the  large  scale  process  is  Gaussian  and  the  small  scale 
is  exponential  gives 
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The  characteristic  function  expansion  is 
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where  the  n=0,  J=0  term  is  again  omitted.  For  Gaussian  large  and 
small  scale  correlation  functions,  Eq.  (25)  is  used,  and  Eq.  (58) 
becomes 
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For  exponential  large  and  small  scale  correlation  functions  Eq.  (28) 
is  used,  and 
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The  incoherent  scattering  coefficients  for  an  exponential  large  scale 
process  and  a Gaussian  small  scale  process  can  be  obtained  from 
Eqs.  (59)  and  (to)  by  interchanging  (g^./^)  and 


Preliminary  numerical  results  were  obtained  from  the  formulas 
given  above  and  were  compared  to  data  taken  at  one  frequency  from  the 
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smoothest  of  ARL's  rough  model  surfaces.  The  data  were  taken  at 
j .00  kHz  from  Surface  I,  which  has  an  rms  height  of  O.091  in.  and 

I correlation  lengths  of  1.4  in.  for  Gaussian  correlation  functions  and 

2.55  in.  for  exponential  correlation  functions.  Figure  1 shows 
results  for  Gaussian  large  and  small  scale  distribution  functions  with 
Gaussian  large  and  small  scale  correlation  functions.  From  the  graph, 
it  is  obvious  that  the  basic  shape  of  the  theoretical  curves  is  correct 
and  that  the  addition  of  the  small  scale  roughness  gives  the  necessary 
flattening  effect  at  low  grazing  angles.  However,  it  is  also  obvious 
that  the  fit  to  the  experimental  data  is  not  good  in  the  region  of 
50  to  80°  grazing  angle.  It  is  believed  that  this  is  in  part  due  to 
the  approximate  treatment  of  the  slopes.  Similar  results  are  seen 
in  Fig.  2,  using  an  exponential  distribution  for  the  large  scale 
roughness,  and  in  Fig.  5,  using  aji  exponential  correlation  function. 

In  addition,  use  of  the  exponential  correlation  causes  the  large 
values  near  normal  incidence  and  somewhat  obscures  the  effects  of  the 
small  scale  roughness. 

The  slight  rise  in  all  of  the  curves  at  low  grazing  angles 
would  be  eliminated  by  including  the  effects  of  shadowing.  Although 
for  some  of  the  cases  shown,  the  fit  of  theory  to  experiment  is  net 
bad,  it  should  be  noted  that  there  are  two  major  factors  which  indicate 
that  further  refinement  of  the  theory  is  necessary.  The  first  of 
these  factors  is  that  fairly  large  values  of  the  small  scale  rms 
height,  a„,  were  necessary  to  produce  the  curves  shown.  The  second 

^ 4 

is  that  the  overall  k frequency  dependence  will  not  allow  fits  to  be 
made  to  the  data  at  other  frequencies  without  significantly  changing 
the  values  of  the  input  statistical  parameters.  It  will  be  shown  that 
both  of  these  problems  can  be  attributed,  in  part,  to  the  slope 
approximations.  In  the  next  section  the  Fresnel  model  will  be  developed 
for  composite  surfaces  using  the  exact  slope  treatment. 
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LARGE  SCALE  ROUGHNESS 


SURFACE  I 
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The  total  scattering  coefficient  for  Gaussian  large  and  small 
scale  distributions  and  Gaussian  correlation  functions  will  now  be  put 
into  a form  that  is  easily  compared  with  the  results  of  the  Fresnel 
model  given  in  the  next  section.  The  total  scattering  coefficient  is 
obtained  by  adding  Eq.  (26)  to  the  coherent  scattering  coefficient 
given  by 
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If  g^»l  and  g^«l,  the  summation  over  j can  be  replaced  by  the  two 
lowest  order  terms,  so  that 
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However,  the  first  summation  is  the  incoherent  part  for  the  lai-ge 
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scale  roughness  alone  with  an  extra  factor  of  e , so  that 
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Assuming  that  I »t^,  restriction  of  the  formula  to  the  backscatter 
^ / 2 "(gT*'g?)\ 

case  where  the  coherence,  (p  e j , is  negligible  gives  the  simple 
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Since  g2=k^7^h2  and,  for  backscatter  7=2sin6r  and  F=l/sin6r,  Eq.  (45) 


becomes 


2 _ . ^^W"'lO^  4^2,2\^‘\y^2A 
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This  equation  will  be  compared  with  the  results  of  the  Fresnel  model 
in  the  next  section  in  order  to  determine  the  effects  of  the  slope 
approximations . 


C.  Modified  Fresnel  Model  with  Exact  Slope  Treatment 


It  has  been  shown  previously  that  the  expression  for  the 
scattered  intensity  developed  from  Eq.  (l)  may  be  expressed  as 


where  the  "O"  subscript  indicates  the  zero  slope  term,  and  the  1 and 
"2"  subscripted  terms  arise  as  a result  of  the  exact  slope  treatment. 
When  a realistic  beam  approximation  and  Fresnel  phase  approximation  is 
made,  the  individual  terms  in  Eq.  (47)  are  given’  by 


ikax  b‘ 


ksinQ  cosS 


ikax  B 


where 


6^  is  the  receiver  grazing  angle 
9^  is  the  source  grazing  angle, 
y = sin9.+sin6  , 


A = Jt  ct  P;  the  elliptical  insonified  area. 


\k 


h 


k is  the  wave  number, 


is  the  surface  to  source  distance, 


is  the  surface  to  receiver  dist;uice, 


K = 3/(^  log^Qo),  a beam  parameter, 


^ " ^oo^ioA^oo'^^10^’ 


^ = "^00^10/ ®i^- 


k'^a^'/aKE^+K/aa^, 
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N = k‘"p^/2KR^+K/2p^, 


and 
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y represents  the  average  intensity  modulation  due  to 
the  surface  roughness . 


For  a Gaussian  distributed  surface  with  two  independent  roughness 
scales,  the  characteristic  function  is  given  by  Eq.  (20) 


'-ik7j^(^^+?2)-(?|+^pj\^  -g^(l-c^)  -ggCl-Cg) 


(51) 


where 


= k“^7 


2 2.  2 


gg  = k 7 hg, 


h^  is  the  rms  value  of  the  large  scale  roughness, 


hg  is  the  rms  value  of  the  small  scale  roughness,  and 


Cg,c^  are  the  normalized  surface  autocorrelation  functions  for 


the  two  scales  of  roughness. 

In  order  to  insure  that  the  results  will  be  simple  we  shall  make  two 
assumptions : 


» 1 


and 


(52) 


gg  « 1 
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With  these  assumptions  we  may  make  the  "large  g"  approximation  in  the 
large  scale  characteristic  fiinction, 

-gl(l-Ci)  -g^l c^(0)| (x^+y^)  (55) 

e = e 

and  the  "small  g"  approximation  in  the  small  scale  characteristic 
function, 


e 


-gp 


(5^) 


To  a reasonable  approximation,  c^  is  given  by  a Gaussian  correlation 
function.  Since  no  physical  knowledge  of  the  small  scale  roughness 
is  available,  it  will  be  arbitrarily  assumed  that  it  is  also 
Gaussian;  thus 
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)/ 1 ^-gg  ■ (®1  l"^  V * 2)^  "®2 

e + g2  e 


(55) 


„ r 2 2./  2 , 2 2J  .2 

-g^(x  +y  )/^.,^-g2  -(x  +y  )/^2"®2 

= e + g2  e , 

In  Eq.  (55)  it  has  been  assumed  means  that  for 

moderately  large  frequencies  the  slopes  of  the  small  scale  roughness 
must  be  much  larger  than  those  of  the  large  scale  roughness. 


Additional  quantities  which  are  required  for  Eqs.  (48)  through 
( 5^'' ' are 
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(56) 


and 


Now  defining  the  scattering  coefficient  in  the  usual  manner, 

0=<l2^^rQQ+r^Q^^,  examinaticn  of  the  first  terms  in  the  brackets  in 

Eqs . (58)^  (59)^  (60)  reveals  that  these  are  the  terms  which  arise 

“®2 

solely  from  the  large  scale  roughness,  except  for  the  factor  e 
The  terms  in  the  second  set  of  brackets  arise  solely  from  the  small 
scale  roughness.  Regrouping  terms,  simplifying,  and  denoting  the  large 
scale  scattering  coefficient  by  , the  scattering  coefficient 
obtained  from  Eqs.  , (59)>  and  (60)  becomes 
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(61) 


The  scattering  coefficient  given  by  Eq.  (61)  contains  several 

interesting  featvires.  First,  the  coefficient  may  be  represented  as 

the  sum  of  the  effects  due  to  the  large  and  small  roughness.  For  most 

“gp 

frequencies  of  interest,  e =1;  thus  the  two  effects  axe  practically 
independent  of  one  another.  Secondly,  the  angular  dependence  of  the 
two  effects  is  markedly  different.  The  scattering  coefficient  for  the 
large  scale  roughness,  starts  out  at  relatively  high  levels  and 
decreases  rapidly  with  decreasing  grazing  angles,  while  the  small 
scale  scattering  effect  is  essentially  constant  at  the  low  grazing 
angles  where  it  becomes  dominant.  Finally,  it  should  be  noted  that 
the  heuristic  model  and  the  modified  Fresnel  model  of  the  scattering 
process  give  identical  results  for  the  small  scale  scattering.  This 
is  because  these  models  differ  only  in  their  treatment  of  the  inco- 
herent component  near  the  specular  direction.  When  the  incoherent 
component  dominates,  as  in  the  case  of  small  angle  backs cattering, 
the  results  of  the  two  models  are  equivalent,  thereby  confirming  the 
suitability  of  the  Fraunhofer  phase  approximation  for  strongly  inco- 
herent fields. 


During  the  next  quarter,  the  literature  will  be  searched  for 
data  on  large  and  small  scale  roughness,  sea  surface  spectra,  and  the 
relation  of  the  resulting  statistics  to  measurable  environment  parame- 
ters such  as  wave  height  and/or  wind  speed.  This  will  allow  model 
input  statistics  to  be  defined  in  terms  of  the  environmental  parameters. 
In  addition,  the  backscattering  model  will  be  modified  to  include  time 
dependence . 
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III.  CONVERGENCE  ZONES 


A.  Introduction 

Convergence  zones  are  regions  of  focused  sound  from  a shallow 

source  in  deep  water.  The  acoustic  waves  are  focused  at  or  near  the 

surface  hy  refraction  at  great  depths.  Convergence  zones  occur 

typically  at  range  intervals  of  80  kyd  or  less  when  the  near  bottom 

sound  velocity  is  greater  than  the  sound  velocity  near  the  surface. 

In  general,  the  convergence  zone  interval  depends  on  the  velocity- 

2 

depth  profile  and  the  location  of  source  and  receiver.  For  example, 
in  the  Atlantic  and  Pacific  Oceans,  the  critical  depth  is  between 
kOOO  and  5000  yd.  The  resulting  convergence  zones  appear  every  60  to 
do  kyd.  When  the  critical  depth  is  around  2000  yd  (as  in  the 
Mediterranian  Sea^)  the  zones  appear  every  30  to  40  kyd. 

Theoretical  treatments  of  convergence  zone  phenomena  are  based 

on  both  ray  and  wave  methods,  but  most  examinations  have  used  ray 

theory  techniques.  Useful  summaries  of  ray  theory  methods  for  studying 

4 5 

convergence  zones  have  been  given  by  Hale  and  Arase  . 

Ray  tracing  yields  a qualitative  picture  of  the  distribution 
of  energy,  and  it  also  provides  a way  of  observing  easily  the  changing 
characteristics  of  convergence  zones  when  the  profile  parameters  and 
source  depths  are  varied.  A typical  ray  trace  is  shown  in  Fig.  4. 

The  source  was  at  a depth  of  100  yd  and  the  velocity-depth  profile  was 
specified  by  the  Epstein  function.  Note  that  the  critical  depth  was 
approximately  2000  yd  and  that  the  convergence  zones  occur  about  every 
52  kyd. 


I 


Ray  methods  can  also  give  range  predictions,  hut  problems  arise 
for  intensity  predictions  due  to  the  caustics  (a  caustic  is  defined 
as  the  locus  of  focal  points  at  which  adjacent  rays  touch  each  other). 
At  caustics,  ray  theory  intensity  calculations  go  to  infinity  because 
(SR/'^)  goes  to  zero.  In  order  to  make  ray  theory  intensity  calcula- 
tions in  caustic  regions,  Pedersen  and  Gordon^  included  various 

diffraction  corrections  in  their  calculations.  They  found  that  the 

7 

correction  derived  by  Brekhovskikh  was  superior  to  the  correction 

8 

derived  by  Marsh.  These  diffraction  corrections  allow  ray  theory 
intensity  calculations  to  be  made  in  both  caustic  regions  and  in 
shadow  zone  regions . 

When  normal  mode  methods  are  applied  to  convergence  zone 
propagation,  there  are,  in  general,  computational  problems  due  to  the 
large  number  of  trapped  modes . There  is  also  the  problem  of  finding 
a velocity-depth  profile  which  permits  the  depth-dependent  differential 
equation  to  be  solved  in  terms  of  functions  which  can  be  calculated. 
Figixre  5 shows  a wave  theory  plot  of  transmission  loss  versus  range 
for  the  Epstein  profile  of  Fig.  U.  The  frequency  was  10  Hz  with 
15  trapped  modes.  For  higher  frequencies,  many  more  modes  are  trapped 
and  the  usual  computational  problems  arise.  There  is  however  qualita- 
tive agreement  between  Figs . U and  5 • 

Since  ray  theory  and  conventional  wave  theory  are  both  limited, 

several  authors  have  attempted  using  WKB  (Wentzel-Kramers-Brillouin) 

9 

techniques . Bartberger  and  Ackler  have  been  able  to  use  WKB  methods 
for  convergence  zone  calculations  at  frequencies  up  to  JOO  Hz . Also 
Leibiger  and  Lee^^  have  applied  WKB  methods  for  convergence  zone 
calculations  up  to  5 kHz . 

In  this  section  a brief  outline  will  be  given  of  Leibiger 's 
method  utilizing  a more  general  Green's  function.  The  basic  idea  is 
to  transform  the  usual  slowly  converging  normal  mode  expansion  into 
an  alternate,  more  rapidly  converging  series.  Such  a representation 
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way  derived  by  Brernmer  by  considering  the  total  field  to  be  made  up 
r)f  a sum  of  multipaths  or  hops.  The  ray  or  hop  series  was  obtained  by 
expanding  the  integrand  (the  inverse  of  the  Wronskian)  of  the  basic 
contour  integral.  Individual  integrals  were  then  evaluated  to  yield 
the  hop  series.  Leibiger  rederived  Brernmer 's  work  in  his  doctoral 
dissertation  and  has  applied  his  techniques  to  many  of  the  propagation 
modes  encountered  in  the  ocean.  In  the  next  section  the  Bremmer-Leibiger 
method  will  be  discussed  using  Green's  function. 

B.  Green's  Function  Expansion 

12 

The  basic  contour  integral  is  given  by 

P(r,z,z^)  = y G(z,z^)  H^(|r)|d§  , (62) 

where  P represents  the  acoustic  pressure,  r is  the  range  coordinate, 
and  z is  the  depth  coordinate.  The  Green's  function  G(z,z^)  is  given 


’n  (z  )n  (z  )■ 
' n (z  ) 


G(z,z^)  = 


n (z  )n^(z  ) 

^2(z>)+V^2)  - nTTzg-]—- 


o)(n2,n^)  1-R^(z^)R2(z2)  n^(zg)nga^' 


where  the  reflection  coefficients  R^(z^)  and  are  defined  at  the 

boundaries  z^  and  z^.  Other  fiinctions  and  parameters  are  defined  in 
Ref.  12  and  are  omitted  for  brevity. 
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When  the 
Eq . (65) , the 


W?vB  approximation  ia  applied  to  the  denominator  of 
ratio  n^(z^)n^(zg)/n^(zy)n2(zj^)  is  replaced  by 


{6h) 


The  relevant  function  to  be  expanded  in  Eq.  (63)  thus  becomes 


-1 

(65) 

When  the  expansion  of  Eq.  (65)  is  used  in  the  fundamental  contour 
integral  of  Eq.  (62),  the  acoustic  field  is  represented  by  an  infinite 
sum  of  integrals . In  and  near  any  particular  convergence  zone  only 
a small  subset  of  these  integrals  contributes  in  a significant  way  to 
the  total  sum.  A great  saving  of  effort  is  achieved  by  evaluating  the 
integrals  rather  than  the  mode  sum,  especially  at  intermediate  and  high 
frequencies . The  integrals  are  so  ordered  that  each  contributing 
integral  can  be  associated  with  a particular  multipath  of  the  con- 
vergence zone. 

Next  quarter  the  multipath  integrals  defined  by  the  series  of 
Eq.  (65)  will  be  given  explicitly  and  discussed. 
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I’  ABSTRACT 


This  quarter  the  reflection  and  scattering  models  were  further  developed, 

and  a review  of  the  literature  preliminary  to  the  development  of  a convergence  zone 

propagation  model  was  completed.  (U-FOUO) 

In  Section  II,  expansion  of  the  heuristic  scattering  model  to  include  composite 
(both  large  and  small  scale  roughness) surface  statistics  is  given  which  sho'uld 
improve  bac’xscattering  predictive  capabilities.  This  model  does  not  include  slope 
correction  terms  but  has  the  advantage  of  simplicity.  The  more  rigorous  theoretic 
model  which  includes  slope  corrections  terms  is  also  expanded  to  allow  for  composi 
surface  roughness.  (U-FOUO) 

A survey  of  the  literature  covering  convergence  zone  propagation  is  suT.marized  in 
Gection  III.  An  approach  is  outlined  for  the  development  of  a convergence  zone 
".odel  using  a Green's  function  expansion.  This  technique  results  in  an  ordered 
integral  expansion  where  each  integral  can  be  associated  with  a particular  multi- 
path. (U-FOUO ) 


DD.TvM473  > 

S/N  0101.807.6801 


iii^ 


'TFTrn 


curitv^davtiflCBUon 


I 


